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Abstract 

We study isovector collective excitations in nuclear matter by employing the 
linearized Landau- Vlasov equation with and without a non-Markovian binary 
collision term at finite temperature. We calculate the giant dipole resonance 
(GDR) strength function for finite nuclei using Steinwedel-Jensen model and 
also by Thomas-Fermi approximation, and we compare them for "^"^^Sn and 
^^^Pb with experimental results. 
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Giant resonances, in particular giant dipole resonance (GDR), built on highly excited 
nuclear states have been the subject of many experimental and theoretical studies in recent 
years [1~3]. The damping properties and the excitation energy dependence of GDR width 
are still among the open problems in nuclear collective dynamics. There are essentially two 
different theoretical approaches to this problem. The first one explains the temperature 
dependence of the width by the coherent mechanism due to adiabatic coupling of the col- 
lective state with thermal surface deformations [4]. In the second approach, referred to as 
coUisional damping, the coupling with incoherent two particle-two hole states plus a Landau 
damping form the mechanism for temperature dependence [5-11]. 

Investigation of the GDR strength function has been carried out in the mean-field ap- 
proximation without the coUisional damping in [10]. In a previous work, we investigated 
the GDR strength function in infinite nuclear matter in quantal framework by employing 
the linearized version of the extended time dependent Hartree-Fock theory (TDHF) with a 
non-Markovian binary collision term [12]. We then applied our results to finite nuclei using 
the Steinwedel- Jensen model. However, this model is not very reliable for the treatment 
of collective dipole oscillations of finite nuclei since in this model it is assumed that the 
nuclear surface remains constant and the density oscillations obeys a wave equation with 
the boundary condition that the radial velocity vanishes on the spherical surface with ra- 
dius R — RqA^/^. On the other hand, semi-classical approaches based on Thomas- Fermi 
approximation which are easier to handle than quantal calculations have been very useful 
to tackle the problems related to finite nuclei [13]. Indeed, the long wave length limit of 
TDHF equation is the Landau- Vlasov equation which is valid if the spatial variations are 
slow. Moreover, if Thomas-Fermi method is used to determine the nuclear density then the 
Landau- Vlasov equation forms a better approximation to quantal theory even if the spatial 
variations are not small [13]. 

In this work, we study the GDR strength function of finite nuclei by employing the 
Landau- Vlasov equation. We do not consider the coherent mechanism and investigate the 
temperature dependence due to coupling of the collective state with incoherent two particle- 
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two hole states using a non Markovian collision term in the linearized Landau- Vlasov equa- 
tion. We apply our results to finite nuclei using Steinwedel- Jensen model. We then consider 
giant dipole excitations in finite nuclei by expanding our results using Thomas- Fermi approx- 
imation. This way, we not only asses the effects of the collision term but we also compare 
the differences between infinite nuclear matter and finite nuclei results. 

The transport equation [14], for the particle phase-space density f{f,p,t) with collision 
term is the Boltzmann equation 

^f{f,p,t)+Vp e{f,p,t).Vrf{f,P,t) - Vr e{f,p,t).Vpf{f,p,t) = K{f) . (1) 

In the framework of Fermi liquid theory defining the quasiparticle velocity asv — Vp e(f , p, t) 
and assuming that the potential energy U in the Hartree-Fock Hamiltonian is local, one 
immediately obtains the transport equation with collision term for proton [q = p) and 
neutron {q — n) distribution functions fq 

p,-{K + Vq,fq}r,,^^K{fq) (2) 

or 

+ Vr^rf, - V. + y,).Vp/, = K{fq) , (3) 

where h — T + U is the mean-field Hamiltonian, T is the kinetic energy, V is the external 
field, and K(f) is the non-Markovian collision term. When / and U change by a small 
amount around the equilibrium, we then have 

fir,p, t) = /e,(e,) + df{f,p, t) , Uif, t) = Uo{r) + dUif, t) (4) 

with 

1 fdU\ p^ 

The equation of motion of the small amplitude vibrations in the semi-classical limit for 
infinite nuclear matter is then obtained as 
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^5/ + v.VrSf - Vr [SU + 2SV].Vpfe,{ep) = 5K (6) 

which is the well-known the linearized Landau- Vlasov equation with a collision term. In 
this equation 5f = 5fn — Sfp , 6U = 6Un — 6Up and 6V = SVn — SVp {SVn,p = ±dV) are 
differences between the indicated neutron and proton functions. 
The isovector mean field 5U{r,t) can be expressed as [7,11] 

5U{r, t) = /o5p(f, t) (7) 

where /o = Fq{T)/N{T) is the quasiparticle zero-order interaction amplitude, Fq{T) is the 
isovector Landau parameter 



Fi(T) ~ Fi(T = 0) 



12 KeFJ 



(8) 



is the density distribuiton function, = 2 is the spin degeneracy factor and 

is the thermally averaged density of states. For T = A^(0) is given as A^(0) — gpFm/2TT'^h^, 
with pf Fermi momentum. 

Now, we present the details of the calculation of GDR response function for infinite and 
finite nuclear matter by using the linearized Landau- Vlasov equation without and with a 
non-Markovian collision term. The solution of Eq. (11) can be found in form of a plane 
wave for infinite nuclear matter 

5fir,p,t) = Sh,Up) e^[^^-^-(-+^^)*l , (11) 

where r) is the vanishingly small positive number corresponding to an adiabatic switching of 
the field at time t = —oo. From the coUisionless Landau- Vlasov equation, we have 

— * 

Sf{r,p,t) + J^''^^''\ SU + 25V]^0 , (12) 

Oi) + IT] — K.V oep 
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and then integrating / j^^s with weight 1, we obtain 

Sp+[foSp + 2SV]x^'\k,u;)^ 



(13) 



where 



gdp k.v dfeq{ep) 



(14) 



cu + irj - k.v dcp 

is the unperturbed Lindhard function. We can then write response of the colhsionless system 
for an external field 5V oc (^\k.f-{u+iv)t] ^ 

Sp 2x^^\k,u) 



n^{k,uj) = — 



(15) 



Performing integrations in Eq.(14), we can find the real and imaginary parts of x^^\^-i'^) ^ 
(for details please refer to [7,11]) 



Imx^^\kM--\N{Q)s 



kpF 



(16) 



and 



Rex^'\k,u) 



mw 



kpF 



/eg(s^e) 



(17) 



Here, 



3peg J (27rfi) 



_gdp_ 



(27r^)3- 

are quasiparticle average kinetic energy and density, respectively and 



(18) 



(19) 



mix) ep 



1/2 



(20) 



kpF V e 

The strength distribution function is obtained from the imaginary part of the response 
function [15] 

2/mx^^^(A;,cu) 



5(A^,w) = --7mn°(A^,u;) = 

^ ^{\-rhReX^%Ku)Y-r{hImx^^\Kuj)Y 



(21) 



The solution of the hnearized Landau- Vlasov equation with coUision term for the infinite 
nuclear matter is given as 

-i(u; + ir])5f + ik.v5f -ik.v^^^^[foSp+2SV]^5K , (22) 

from which we obtain 

5p+[foSp + 25V]x^'\k,uj) = -[fo5p + 25V]x^'\k,uj) . (23) 

The coUisional response function x'^'^\k,(jj) can be expressed as [12] 

(27rAi)'' J \ 2 J TT w — £3 — £4 + ei + £2 + ^r7 



where fi = feq{ei), f^^l- feq{ei), ^ Qi + Q2 - Qs - Qa with Qi = 1/ 
ej = {m/2)vf and W(12;34) denotes the basic two-body transition rate 



w — k.Ui 



^^'2; 34) = < ^kl^ > mp.+P2-Ps-p.) (25) 

which can be expressed in terms of the scattering cross-section as 

Then the retarded response function with collision term is obtained as 

U-rk,-) = -^,-^^^^^ (27) 
SV l + f,x{k,u;) 

with xik,uj) = x*^^''(A;, a;) + x^'^Hk,uj). Thus we can rewrite the strength distribution for our 
system with collision term as 

sCk,w) = -'-ImU-^%u) = -i 2ImxCk,u^) ^ _ 

^ ^{l + foRex{k,u;)y + {foImx{k,u;)y 

The strength function satisfies the following energy weighted sum rule [10] 

djjujSik^uj) — - — po (29) 
2m 

where po — 0-lGfm~^ is the saturation density of infinite nuclear matter. 
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In the calculations of the coUisional response function, we use conservation laws and 

symmetry properties. It is possible to reduce the twelve dimensional integrals to five fold 

integrals by the total and relative momentum transformations (P = pi +p2, P' = Ps + Pa, 

and relative momenta q — {pi — p2)/2, q' — {ps — p4)/2) before and after the collisions. We 

neglect the real part Rex^'^\k, u) of the function x^^^ {k, u) in our calculations. The coUisional 

response function x^^^ (k, cu) has a singular behavior arising from the pole of the distortion 
r ~* "I 

functions, Qi — 1/ w — k.vi . We avoid this singular behavior by incorporating a pole 
approximation. In the distortion functions, we make the replacement a;^a;c — ir/2 where 
ujd and T are determined from 1 + /ox*^^^ (A;, a;) = at each temperature that is considered. 
So, we evaluate the remaining five dimensional integrals numerically by employing a fast 
algorithm. In the evaluation of momentum integrals we make the replacement {da/dfl)^^ —>■ 
c"pn/47r with = 40 mb, thus neglecting the angular anisotropy of the cross section. 

In order to apply our results to nuclear dipole vibrations and finite nuclei, we work within 
the framework of Steinwedel and Jensen model which describes the GDR in heavy nuclei as 
a volume polarization mode conserving the total density Po = Pn + Pp ior infinite nuclear 
matter [16] where neutron and proton oscillate inside a sphere of radius R as 

Pp{f,t)-pn{f,t)<xsm{k-r)e'''' . (30) 

According to this model, we choose the wave number of the normal mode as k — Tr/2R. We 
apply Steinwedel and Jensen model to GDR in ^'^^Sn and ^^^Pb, and we take R = 5.6 fm 
k = 0.28 fm-^ for ^^o^^ and P = 6.7 fm A; = 0.23 fm'^ for ^^^Pb according to R ^ 1.13^^3. 

The Landau parameter Fq(T — 0) can be expressed as a function of the symmetry energy 
coefficient a-,- in the Weizscker mass formula at zero temperature as follows [17] 

F^(T = 0) = ^-1 . (31) 

For the value of = 28 MeV we have F^{T = 0) = 1.33. The value of the F^{T) decreases 
with temperature because of the decrease of the thermally averaged level density N(T) [7]. 

So far, our GDR calculations have been for infinite nuclear matter. In the rest of the 
paper we will employ the Thomas-Fermi approximation (TF) to calculate GDR response 



function by using the linearized Landau- Vlasov equation with and without the coUision 
term at finite temperature for finite nuclei. The Thomas-Fermi theory, together with its 
extensions, is the semiclassical treatment of nuclear dynamics in its independent particle or 
Hartree-Fock approximation and can be explained from quite different points of view [15,16]. 
We evaluate the expression for the GDR in the TF approximation, which corresponds to a 
semi-classical transport description of the collective vibrations. The solution of the linearized 
Landau- Vlasov equation with and without collision term 

^5/ + vVrSf - V, [SU + 2SV].^Mep, r) = SK (32) 

for finite nuclei is obtained with the local plane wave ansatz 

Sf{f,p,t) = 6h,^{f,p) e^[^^--(-+^'')*l . (33) 

Response of the finite system without and with collision term is then obtained as 



where XTFik,uj) = Xtf{^)^) + Xtf(^)^)- The strength distribution function without and 
with collision term are 

STF{k,w) = --Imn%p{k,uj) , STF{k,w) = --ImW^^Xk.uj) (35) 

TT TT 

with 

Am dfpir) x^\k,^,r) (36) 

where i — 1,2 , and 

fr = \jdrp{r)Ur) . (37) 

The function x^"^^ {k, oj,r) is obtained by evaluating the collision term given in Eq. (24) using 
Thomas-Fermi approximation. We determine the nuclear density p{r) for the finite nuclear 
matter in TF approximation using a Wood-Saxon potential with a depth Vq = 44 MeV , 
thickness parameter tp — 0.67 fm and sharp radius R — 1.13^4^/^ [16], 



V{r) = ^ , p{r) = —kF{r)'Q{X-V{r)) , k^ir) = (-^[V{r,) - Vir)]) 

(38) 

where Tc is the critical radius for a mass number A is defined as 

r drp(r) . (39) 
Jo 

Here A = ^(^c) a-nd, the expression for can be numerically integrated for a given 

A to determine Tc- For finite nuclear matter, the interaction amplitude is fo{r) = 
3V^'(r)p(r)/(2eF(r) N{r,T = 0)) which is related to the parameters of the simplified Skyrme 
force as 

Vi{T)^-\h{x^ + \)-\hp{T) . (40) 

We use the following parameters: — —983.4 MeVxfm^, ^3 = 13106 MeVxfmP, and 
xo = 0.48 [10]. 

We show our results for the GDR strength function with and without the collision term in 
Fig. 1 and in Fig. 2 for and ^°^P6, respectively, calculated using infinite nuclear matter 
formalism within the framework of Steinwedel- Jensen model. In Fig. 3 and in Fig. 4 we show 
the GDR strength function with and without collision term for ^^°5'n and ^°^P6, respectively, 
calculated employing Thomas-Fermi approximation. In this figures, we also compare our 
results with the normalized experimental data taken from [1]. The temperature parameter T 
in the mean occupation number functions /(e, T) is related to the experimental temperature 
T* asT — T* ^ dp where ap — Air^/Aep is the Fermi gas level density parameter and ap 
is the energy dependent empirical level density parameter [1] . 

From these figures, we first note that without the collision term the position of the peak 
of the strength function does not change appreciably with temperature. This behavior is 
in accordance with the experimental results [1,2]. Since we neglect the real part of the 
coUisional response, when we include the coUisional term the average position of the peak 
values of the strength functions do not change again, but as the result of the collisions the 
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overall shape of the strength function changes somewhat and this change becomes more 
pronounced with increasing temperature. In the case of infinite nuclear matter this change 
has the tendency to improve the agreement with experimental results for ^"^^Sn but this 
tendency is much less pronounced for ^°^P6 as it can be seen by comparing Fig. 1 with Fig. 
2. On the other hand, for this case of finite nuclei calculations employing Thomas- Fermi 
approximation the change produced by the addition of the coUison term can be clearly noted 
in Fig. 3 for ^^^Sn and in Fig. 4 for '^^^Ph. Indeed, in both cases, the overall agreement with 
the experimental results is much more improved when the binary collision term is included. 

In our work, we obtain a reasonable description of the giant dipole excitations in ^"^^Sn 
and ^°^P6 using semi-classical approach with Thomas-Fermi approximation as compared to 
infinite nuclear matter formalism, and we demonstrate the importance of the collision term 
which improves the agreement of the calculated strength functions with the experimental 
results. We believe that inclusion of the coherent damping mechanism into our formalism 
will extend our description further. 
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FIG. 1. The GDR strength function of ^^°Sn obtained using Steinwedel- Jensen model. Sohd 
and dashed lines show the response function without and with the collision term, respectively. The 
normalized data is taken from [1]. 
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FIG. 2. The GDR strength function of ^'^^Pb obtained using Steinwedel-Jensen model. Solid 
and dashed lines show the response function without and with the collision term, respectively. The 
normalized data is taken from [1]. 
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FIG. 3. The GDR strength function of ^^°Sn calculated by Thomas-Fermi approximation. Solid 
and dashed lines show the response function without and with the collision term for the finite 
nuclear matter. The normalized data is taken from [1]. 
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FIG. 4. The GDR strength function of ^°^Pb calculated by Thomas-Fermi approximation. 
Solid and dashed lines show the response function without and with the collision term for the finite 
nuclear matter. The normalized data is taken from [1]. 
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